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Abstract. Radiation-hydrodynamical simulations of surface convection in low-mass stars can 
, be exploited to derive estimates of i) the efficiency of the convective energy transport in the 
• stellar surface layers; ii) the convection-related photometric micro-variability. We comment on 
the universality of the mixing-length parameter, and point out potential pitfalls in the process 
of its calibration which may be in part responsible for the contradictory findings about its 
variability across the Hertzsprung-Russell digramme. We further comment on the modelling of 
, the photometric micro- variability in HD 49933 - one of the first main COROT targets. 
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d . 1. Introduction 

Radiation-hydrodynamical (RHD) simulations of surface convection in low-mass stars 
have reached a high level of maturity. Such simulations provide quantitative predictions 
^ , about the spatial and temporal statistics of the flows taking place in the stellar surface 
layers. We used the 3D RHD code CO''' BOLD to study convective flows in late-type 
stars at different metalicity. This contribution deals with two distinct applications of 
P\| CO^BOLD models. The first one is related to the efficiency of the convective energy 
transport in convective envelopes, the second one to the low-level photometric variability 
' related temporal evolution of the surface granulation pattern. We shall rather discuss 
(yO , problems than solutions with a focus on dwarfs. 



2. Convective energy transport in the envelopes of late type-dwarfs 

It is well known from the theory of stellar structure that convection is generally an 
' efficient means of transporting energy, and it establishes a thermal structure close to adi- 
abatic. Only in vicinity of the boundaries of convective regions noticeable deviations from 
adiabaticity occur. In convective envelopes of late-type stars the upper boundary of the 
convective envelope - usually located close to or even in the optically thin layers - con- 
stitutes the bottle-neck for the energy transport through the stellar envelope assigning a 
special role to it. Despite it is geometrically thin and contains little mass it largely deter- 
mines the properties of the convective envelope as a whole. It is the value of the entropy 
of the adiabatically stratified bulk of the convective Sonv which is most important from 
the point of view of stellar structure since it influences the resulting radius and effective 
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temperature of a stellar model. Scnv is controlled by the efRciency of convective and ra- 
diative energy transport in the thin, superadiabatically stratified surface layers. Detailed 
RHD simulations can be applied to model this region allowing to quantify the mutual 
efficiency of the convective and radiative energy transport and predict Senv Comparing 
the RHD predictions to standard ID models based on mixing-length theory (MLT) the 
value of Sonv can be translated into a corresponding mixing-length parameter qmlt (see 
Ludwig et al. 1999 for details). 

In stellar evolution calculations the free mixing-length parameter is usually calibrated 
against the Sun. However, it is unclear whether mixing-length theory provides a suit- 
able scaling of the convective efficiency at constant auur across the Hertzsprung-Russell 
diagram (HRD), and a lot of work was invested to address this issue empirically. Unfor- 
tunately, hitherto, no coherent picture emerged. Here we report on an update of earlier 
work on the theoretical calibration of omlt based on RHD simulations. While not at 
all comprehensive it illustrates some of the pitfalls in the process which may in part be 
responsible of the blurred and sometimes even contradictory picture which emerged so 
far concerning the variability of the mixing- length parameter. 



2.1. Comments on the functional dependence of auLT on stellar parameters 

Discussions which took place during the symposium led the authors to add a comment 
about the question which stellar properties govern the value of the mixing-length param- 
eter aiviLT- We argued above that RHD models of the surface layers are able to provide 
information about Scnv and Q!mlt- The models are characterised by the atmospheric pa- 
rameters and consequently the functional dependence of aMLT can be described in terms 
of them. Whether the standard atmospheric parameters Toff and log g together with the 
chemical composition are the most suitable coordinates is not clear. One might speculate 
that, e.g., the surface opacity is a physically more relevant quantity. Nevertheless, for 
unevolved late-type stars the conditions at the stellar surface govern the global envelope 
structure and the standard atmospheric parameters are suitable coordinates to param- 
eterise them. Global stellar parameters (mass, radius, or age) play merely an implicit 
role. The situation only changes when the size of the granular cells or the thickness of 
the superadiabatic layer become comparable to the stellar radius which might happen in 
giants. To our opinion one should usually avoid to express changes of qmlt in terms of, 
e.g., stellar mass or age since it tends to obscure the underlying physics. 



2.2. Q!MLT from six 3D RHD models 

Squares in Figs. [1] and [2] mark the parameter combinations in the Tcfi-log 5-plane of six 
3D RHD models we used for calibrating aMLT- All atmospheres belong to dwarfs, three of 
solar metalicity, three of a metalicity of 1/100 solar. Our new 3D results are superimposed 
on earlier results obtained from a grid of 2D RHD models by Ludwig et al. (1999), and 



Freytag et al. (1999) There are obvious differences of up to w 0.8 between the 3D and 
2D calibrated values. Does this indicate that a theoretical calibration with the help of 
RHD models is hopelessly inaccurate? 

The first question to answer is what is a big and what is a small difference when it comes 
to Q!MLT- At present interferometry and analysis of eclipsing binaries provide stellar radii 
to an accuracy of about 1 %. Looking at stellar structure models (see, e.g., Lebreton et al. 
2001, their Fig. 10) for late-type main-sequence stars the highest sensitivity of the stellar 
radius and effective temperature to changes of Q!mlt is found at about one solar mass. 



Christensen-Dalsgaard (1997) finds for the Sun a sensitivity of (51n_R « —0.245 In aMLT- 



Hence, considering the present observationally achievable accuracies we would like to 
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Figure 1. Theoretically calibrated mixing- length parameters at solar metalicity in the 
Tcff-log g-plane. Squares mark our new results based on 3D models, the other symbols mark 
earlier 2D results. The obtained values are given by the numbers, the isolines represent a smooth 
fit to the 2D data. 
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Figure 2. Same as Fig. [l]but for 1/100 of the solar metalicity. 



know QiMLT to better than 4%, and the differences between the 2D and 3D results are 
clearly relevant. 

We did not mention yet the methodological changes we introduced when calibrating our 
3D-based omlt's in comparison to the earlier works: i) Our new values were calculated 
with a different mixing- length dialect. The earlier results assumed the formulation given 
by 'Bohm-Vitense (1958) while we now used the formulation given by Mihalas (1978) ii) 
Following the general trend in stellar evolution theory, our ID comparison models are now 
full-fledged stellar atmosphere models instead of integrating a prescribed T(r )-relation to 
describe the atmospheric temperature run as was done in the earlier results. As we shall 
see in a moment this is a crucial point, in particular for the cool, metal-poor model at 
5000 K. iii) It is now well established that the convective energy transport operates more 
efficiently in 3D than in restricting 2D symmetry. Hence, we would expect a systematic 
bias towards higher mixing- length parameters in 3D relative to 2D. 

Having in mind that we expect larger ckMLT-values in 3D we think that the absolute 
differences are perhaps not surprising and now rather focus on differential trends with 
Toff- The statistical uncertainties in the 2D calibration amount to ±0.05. Taking together 
the different flavours of MLT, and the different ID comparison models used in the 2D 
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Figure 3. Entropy profiles for ID stellar atmosphere and stellar structure models for various 
mixing-length parameters in comparison to the mean 3D stratification for the metal-poor model 
at Teff^SOOOK, logg=4.5, and [M/H]=-2.0. All ID models employ the MLT formulation by 



Mihalas (1978) The dashed line depicts the entropy Scnv predicted by the 3D model for the 
adiabatically stratified part of the convective envelope. 



and 3D calibrations the trends start to look similar in 2D and 3D with the exception of 
the coolest metal-poor model which is far off. What is the reason? 

Figure [3] shows a comparison between entropy profiles of ID MLT based models and of 
the horizontally and time- wised (on surfaces of equal optical Rosseland depth) averaged 
3D model. This kind of comparison provides the calibrated aMLT- It is obvious that the 
result depends on the employed ID models. Using stellar atmosphere models produces a 
mixing- length parameter which is almost by 0.5 larger than one based on models which 
use a prescribed r(T)-relation in the convectively stable part of the stellar atmosphere. 
The reason is to a large extend the difference in the entropy minimum Smin attained in 
the deep atmospheric layer in combination with the small overall entropy jump - the 
difference Senv — Smin - in the star which amounts to only « 1/6 of the solar value. Stars 
of higher Tes show larger entropy jumps which makes the exact level of the entropy 
minimum less critical. The quite different differential behaviour of ckmlt in 3D relative 
to 2D is mainly a result of the different choice of comparison model. 

While one might consider this as mere technicality we rather believe that part of the 
confusion about trends and absolute values of the mixing-length parameter have their 
origin - besides observational problems - in different and not clearly specified procedures 
of how the MLT-related quantities are computed and the atmospheric structure integra- 
tion are actually performed in models. It is, e.g., still quite common to use simple grey 
r(T)-relations for describing the atmospheric temperature run. Figure [3] illustrates that 
the particular choice might generate mismatches corresponding to substantial changes in 
the value of the mixing-length parameter necessary to restore the actual scaling of the 
envelope entropy with changing stellar effective temperature. 

Finally, Fig. [3] illustrates that even full-fledged model atmospheres are not always able 
to reproduce the entropy minimum in 3D models closely. The reason is that for metal- 
poor dwarfs 3D models predict large deviations from radiative equilibrium conditions in 
the formally convectively stable part of the atmosphere (in the figure apparent by the low 
entropy of the 3D model at low optical depth). In ID stellar atmospheresradiative equi- 



librium conditions are assumed to fix the atmospheric temperatures. Trampedach (2007) 



puts forward the idea to extract the mean T(r )-relation from 3D models for performing 
the atmospheric structure integration. This obviously provides a better match to the 3D 
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Figure 4. Temperature-pressure profiles of a 3D red giant model (tliick solid line) in comparison 
to ID stellar atmosphere models of different q:mlt leaving out (thin solid lines) or including 
(dashed and dashed-dotted line) turbulent pressure (for details see text). 

atmospheric entropy minimum, and hence a "cleaner" calibration of the mixing-length 
parameter, however, also is more difhcult to implement in existing stellar evolution codes. 
The 3D calibrated values of ajviLT by Trampedach (for solar metalicity and MLT formu- 
lation by Bohm-Vitense) are somewhat lower than our values presented here but show a 
similar trend with Toff on the main-sequence. 

2.3. Turbulent pressure trouble 

In the main-sequence models discussed above turbulent pressure plays generally only a 
small role but becomes relatively more important towards lower gravities - and causes 
extra trouble when one is interested in a well-defined calibration of the mixing-length 
parameter. Figure |4] shows the average temperature profile of a 3D red giant model 
(Toff=3600K, log5=I.O, [M/H]=0.0) in comparison to standard ID model atmospheres 
of the same atmospheric parameters. While turbulent pressure Ptmb is naturally included 
in the RHD simulation it is modelled in a ad-hoc fashion in ID models assuming a 
paramterisation Pturb = /turbP^^c' where /turb is a free parameter of order unity, p the 
mass density and the convective velocity according MLT. 

Figure S] shows that in the 3D models the turbulent pressure "lifts" the temperature- 
gas pressure profile towards lower pressures which is essentially impossible to reproduce 
in the ID models ~ irrespective of the choice of q;mlt and /turb- The failure is related to 
the local nature of MLT confining the action of the turbulent pressure gradients strictly to 
the convectively unstable regions. While formally one can still match the thermal profile 
of the 3D model in the deeper layers by ID profile with suitably chosen amui and/or 
/turb such a match becomes physically little motivated and is unlikely to provide a robust 
scaling with changing atmospheric parameters. An improved ID convection description 
including non-local effects like overshooting is clearly desirable to handle this situation. 
An empirical calibration of q;mlt using giants is likely to suffer from ambiguities related 
to the way turbulent pressure is treated in ID models. 




3. Modelling the photometric micro- variability in HD 49933 

The stochastically changing granular flow pattern on the surface of late-type stars 
causes low-level residual brightness fluctuations in stellar disk-integrated light. Ludwig (2006) 
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Figure 5. Predictions for the temporal power spectral density of the photometric brightness 
fluctuations exhibited by HD 49933. 3D RHD simulations are depicted by solid lines and labelled 
by Toff/log (;/[M/H]. Dashed lines depict the result when adding the estimated photometric 
noise level. For further comparison a spectrum obtained with the COROT light curve simulator 
(Baudin et al. 2007) is shown. 

described a method to obtain predictions of the observable temporal power spectrum 
from local-box RHD simulations if - besides the atmospheric parameters - the radius of 
the target star is known. One of the first primary asteroseismic targets of the COROT 
satellite mission was the metal-depleted (-0.3 ^ [M/H] ^ -0.4) F-dwarf HD 49933 for 
which a high-precision photometric light-curve was obtained. The data is not released yet 
so that we are only showing theoretical predictions obtained from two RHD simulation 
runs at Tcff=6750K, log (7=4.25 and two metalicities [M/H]=0.0 and -1.0 bracketing the 
observed metalicity of HD 49933 and having atmospheric parameters close to the spec- 
troscopically measured values. Figure [5] depicts the power spectra obtained from the two 
RHD simulations. For smoothing, the raw RHD spectra were fitted by a simple analytical 
model, and oscillatory peaks from the acoustic eigenmodes of the computational box were 
removed. Hence, the pure convection-related signal is shown. The observed signal should 
fall between the two model predictions, at least at around 1 mHz where the photon noise 
and the signal due to magnetic activity are expected not to dominate. Unfortunately, at 
present it looks that this is not quite the case. An actual comparison will be presented 
in an upcoming paper by Ludwig & Samadi, so stay tuned. 
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